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Abstract 



We consider the problem of constructing an action functional for physical systems whose classical equa- 
tions of motion cannot be directly identified with Euler-Lagrange equations for an action principle. Two 
, ways of action principle construction are presented. From simple consideration, we derive necessary and 

Oh' sufficient conditions for the existence of a multiplier matrix which can endow a prescribed set of second- 

order differential equations with the structure of Euler-Lagrange equations. An explicit form of the action is 
constructed in case if such a multiplier exists. If a given set of differential equations cannot be derived from 
an action principle, one can reformulate such a set in an equivalent first-order form which can always be 
treated as the Euler-Lagrange equations of a certain action. We construct such an action explicitly. There 
exists an ambiguity (not reduced to a total time derivative) in associating a Lagrange function with a given 
set of equations. We present a complete description of this ambiguity. The general procedure is illustrated 
by several examples. 
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! 1 Introduction 



The problem of constructing an action functional for a given set of differential equations is known in literature 
as the inverse problem of the calculus of variations for the Newtonian mechanics. In its classical setting [1] 
' consists of solving the variational equation 

o: ssM^ 

^ , where gi{t, q^,ct, ...) =0 is some given system of differential equations with respect to unknown functions q*(t), 

^ ■ and S[q\ is a local functional to be determined. The condition of locality requires the existence of a function 

5^ , L{t,q,q...) (Lagrangian) , such that the functional S[q] (action) would be written as an integral 

S[q]^ J dtL . (2) 

In other words, the essence of the inverse problem of the calculus of variations consists of finding a variational 
principle for a given system of differential equations. This problem has been under consideration for more than 
a hundred years. As early as 1887 Helmholtz [l] presented a criterion of commutativity for second variational 
derivatives from which immediately follows the necessary (and with some restrictions also sufficient) condition 
of solvability of the equation ([1]) : 



Sgijt) ^ %(g) 
dqi{s) Sqi{t) 



(3) 



If this condition holds, the system gi{t, q^ ,q\ ...) = is called Lagrangian system, if not the system is non- 
Lagrangian. In 1894 Darboux |2] solved the problem for the one dimensional case. In 1941 the case of two degrees 
of freedom was investigated by Douglas [3]; in particular, he presented examples of second-order equations which 
cannot be obtained from the variational principle. Afterwards many authors (see e.g., 0-[l5] and references 
therein) investigated this problem for multidimensional systems. 
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In the present work we consider the question of the construction of an action principle for a given system 
of differential equations using the integrating multiplier method [3]- [9]. The integrating multiplier is a non- 
singular matrix which being multiplied by a given set of differential equations reduces this set to a standard 
Euler-Lagrange form. In section 2 we present a simple derivation for necessary and sufficient conditions for 
an integrating multiplier for a system of second-order equations. We also construct the expHcit form of the 
Lagrangian in case an integrating multiplier exists. Then we apply our method for investigating the inverse 
problem of some simple models. In particular, we construct an action principle for multidimensional dissipative 
systems. We also consider an example of a linear dynamical system whose equations of motion does not admit 
an integrating multiplier, and, as a consequence, cannot be obtained from the minimum action principle. 

Note that it is always possible to reduce non-Lagrangian second-order equations of motion to an equivalent set 
of first-order differential equations. From the Helmholtz criterion ([3]) we find necessary and sufficient conditions 
for the existence of an integrating multiplier for such equations. It turns out that in the first-order formalism 
an integrating multiplier always exists and can be constructed by means of the solution of the Cauchy problem 
for the equations in question, and this is presented in section 3 and is partially based on results of the works [10] 
and |T5]. Then we construct the action functional explicitly. Thus, we show that systems traditionally called as 
non-Lagrangian are, in fact, equivalent to some first-order Lagrangian systems. As an example, we construct a 
first-order action functional for any linear dynamical system. 



2 Action functional for a set of second-order equations 
2.1 General consideration 

Let a system with n degrees of freedom be described by a set of n second-order differential equations of motion, 
solvable with respect to second-order time derivatives. Suppose such a set has the form 

q'~r{t,q,q)^0, i = l...n, (4) 

where f''(t,q,q) are some functions of the indicated arguments, and by dots above we denote time derivatives 
of the coordinates. Let us construct an action principle for this set. If (0]) cannot be directly identified with 
Euler-Lagrange equations, one can find an integrating multiplier, i.e., a nonsingular matrix hij{t,q,q) that 
being multiplied by |[4)) 

h^,[q^ ~r{t,q,q)]=0 (5) 
reduces this set to standard Euler-Lagrange form for some Lagrangian L{t,q,q), 

^^^^^q^-^q^=0. (6) 
dq^ dtdq^ dq'-dq^ dq^dq^ 

In order to identify ^ with ^ we need to ensure that 



dq^dq^ 



h^j, (7) 
q^=h,,r. (8) 



dq^ dtd(f dq^dq 

Provided that an integrating multiplier is known, equations {Tj-lIHl) can be interpreted as a system of equations 
for a Lagrange function L. We are going to solve the set of equations ll7l)-|[8l). Its consistency conditions will 
give us all the necessary and sufficient conditions for an integrating multiplier. Assuming that L is a smooth 
function of the indicated arguments, the consistency condition for equation ^ imply that 

h -h f9^ 

- ""^^ ' dq^ ~ dq- ■ ^ ' 

If ([9]) does hold one can solve equation ((?]). To this end, we remind that the general solution of the equation 
df /dq^ = gi, provided the vector gi is a gradient, is 



/(<?) = /' 
Jo 



ds q'gii-sq) +c, 
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where c is a constant. Taking the above fact into account, we obtain for L (we do not consider global problems 
which can arise from non-trivial topology of the configuration space) the following representation: 



where 



L ^ K{t,qA) +k{t,q)q' +lo{t,q) , 



K{t,q,q)^ Jdaq^ j dhq\hij{t,q,hqi) 



(10) 



(11) 



qi=aq 



and ^o(i, ?), h{t^ q) are some functions of the indicated arguments. To find these functions, we use equation |[8]). 
Substituting HO]) into (HI), we get 



dK d^K d^K 



dh dh 



dq^ dq^dt dcfdqi \dq^ dqi 
Differentiating this equation over q^ , we obtain: 

dh dl 
dq^ dq 



dh , dlo 



dt dq' 



where 



d^K d^K dhik .,dh,k d 



fdqi dqldqi ' dt ' " dq^ ' dq'' '^^'^^> 



(12) 



(13) 



(14) 



This equation is a differential equation for U. The consistency conditions of equation lfT3|) imply that, first of 
all, the symmetric part of L^fc is zero, which can be written as 



op _ dp 



= 0. 



(15) 



where 



D 



d 



dt 



P 



d 



dqi dqi 



Using l|15p . one can rewrite (fT4l) as follows 

d'^K d^K _^ (u ^P u 

~ dctdq'^ ~ dqf'dq^ + ' " 2 [ dq'' ' dq^ 

Next, Lik does not depend on the velocities, i.e., dLik/dcl = 0, which yields 

dhu _ dhg _ _d_j^ 

Qqi Qqk Qql ' 

And, finally, the Jacobi identity 

dLik , dLki , dLi 



dq' dq' dq'' 







dA,k dAki dA. 



dq' dq' dq'' 



(16) 



(17) 



(18) 



Provided hij obeys equations (fTSj) . (flTl) and (flSl) . k can be found from equation (fT3| . We remind that the 
general solution for k of equation lfT3|) is given by 



kit,q)= daq'^Lki{t,aq) + 



dip {t,q) 
dqi 



(19) 



where (p(t, q) is an arbitrary function. 

Now from equation lfT2|) we can find Zq; to this end let us rewrite it as follows: 

5^ = "^- 



(20) 
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where 

The consistency conditions of (|20|) imply that, first, rrii does not depend on the velocities, i.e., drm/dq'^ — 0. 
This condition is provided by equation ifTSj) . And second, the vector rrii must be a gradient: 

Taking into account (fTSj) and lfT7|) . one gets from l(22|) the following algebraic condition: 

h,,Bi-hk,Bi=Q, (23) 

where 

■ ^ 1 g/" f^df dp 

' 2 9g™ dq^ dq^ dqJ ' 

If hij obeys (|23l) , then from (|20l) one gets 

1 

Zo(i,9)= jdaq^mk{t,aq) +^^P^+c{t), (24) 


where c (t) is an arbitrary function of time. 

Thus, we have proved the following statement: iff for a given set of second-order ordinary differential 
equations ^ there exists a non-singular matrix hij {t,q,q)that obeys equations (E), ( 05)) . flT]) . (G3i and 
f23\} . then this set can be obtained from the variational principle with Lagrangian flC^) . where the functions 
K {t, q, q) , li (t, q) and Iq (t, q) are defined by II 11]) . and \24\) respectively, and the functions ip (t, q) and 

c (t) are arbitrary functions of the indicated arguments. 

The arbitrariness related to the functions ip {t, q) and c (t) enter the Lagrangian (fTO|) via the total time 
derivative of a function F, 

F = Lp{t,q)+ j c{t)dt. 

Note that an integrating multiplier hij , and as a consequence the Lagrange function L does exist, but however 
not for any set of equations ([4]). In Section 3 we consider an example of a dynamical system which does not 
admit the existence of an integrating multiplier. However, if it exists, it is not unique [Il]-[2l], e.g., if the matrix 
hij is an integrating multiplier for a certain set |(4]), it is easy to see that the matrix hij = c hij, where c 7^ 
is a constant, is an integrating multiplier as well. Therefore, Lagrangian ifTOl) leading to the set of equations 
(0]) is not unique; because for this set there exist as many inequivalent Lagrangians as integrating multipliers. 
Lagrangians corresponding to different integrating multipliers are known as s-equivalent Lagrangians. 

In the one dimensional case q — f(t, q, q) — 0, an integrating multiplier is a non-vanishing function h{t, q, q) 
that obeys the equation 

dh dh d , „, . , ^ 

m^'^d-q^d-q^^'^ = '- ^''^ 

This is a first-order partial differential equation which obviously has a solution for any / and initial condition 
h{t = 0, q, q) = ho {q, q). As we can see, an answer to the question whether there exists a solution of the inverse 
problem of the calculus of variations depends on the number of degrees of freedom n. For n = 1 the answer is 
always positive, and there exist as many inequivalent Lagrangians as functions ho {q, q) of two variables. For 
n > 2 the answer is generally negative. 

2.2 Examples 

In this section we consider the possibility of constructing an action principle for some examples of dynamical 
systems. First of all, let us consider dissipative systems. Suppose we have an ideal system with Lagrangian 

• 2 

Lo^^ + V{q), q^{q'}, l=l...n. (26) 
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Let us consider the case when besides the potential conservative force F' = |^ there exist a friction force 

F}„c = aq\ (27) 

where a is a phenomenological friction coefficient which in general can depend on time. The equations of motion 
for such a system have the form 

dV 

+ (28) 

These equations are non-Lagrangian, but for this set it is possible to find an integrating multiplier. In the 
simplest case, when it does not depend on coordinates and velocities, it has the form 

/i,^. =e-2/"'it/jO^ (29) 

where ft^^ is an arbitrary, symmetric, nonsingular, constant matrix commuting with the matrix Vij — d^V/dq^dq^ . 
Using the statement of the previous section, we obtain the following Lagrangian: 



L = -q'h.jq^ + j q'Kj ds. (30) 



If one sets = 5ij , Lagrangian l(30|) can be rewritten as 

L = e-2/"'i*Lo. (31) 

Note that once the friction coefficient goes to zero, Lagrangian l(3T1) transforms into the initial Lagrangian l(26|) . 

Let us now consider the case when the potential in the initial Lagrangian is linear in velocities. For simplicity 
we consider the two-dimensional case 

Lo = ^ (i' + f + f3{xy - yx)) , (32) 



Let us consider this system in the presence of the dissipative force l(30|) . The equations of motion will have the 
form: 

X — ax ~ I3y , 

(33) 

y — fix + ay . 

As was shown in [12] , this system describes a moving charged particle in a uniform magnetic field with radiation 
friction. In this case, 

^ 2a(3 a2_^2 

and from equation l(23l) one immediately gets 

tr(/iy) = /ill + ;i22 = 0. (34) 

It is then easy to find that the general solution of the equations ^ , (fT5|) , (fTZl) , ifTS]) is defined by an arbitrary 
function (/) (^, 77) and has the form 

h - ( ^ + ^ ^{F-F) \ 

where F = (p (^£,e~^* ,^ — e"''*, = x + iy, ^ — a + if3 and the bar denotes a complex conjugation. 
The simplest real solution can be found if we put (j) — We have 

= f ^ y.) . (36) 

x^ + y^ \ y ~x J 
Using the formulas (fTO| . (flTI) . (fT9|) and l(24|) we find the following Lagrangian: 

L = In (i^ +y'^) +y arctan + ax - f3y . (37) 
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The corresponding Euler-Lagrange equations 



XX + yy XX — yy ^ , , 

are equivalent to the initial ones l(33|). with the exception of the point i = y = 0. Thus, we can see that in this 
case the inverse problem of the calculus of variations is solvable. Unfortunately, neither Lagrangian (|37| nor 
any other Lagrangian constructed by the matrix (|35l) in the limit a ^ transforms into the initial Lagrangian 
(|32l) . modulo a total time derivative. This is because, according to the algebraic condition l(34|) . the trace of the 
Hessian matrix of any Lagrangian for the set of equations l(33|) must be equal to zero and this property holds 
true after the limit a — > is taken. On the other hand, the trace of the Hessian matrix of the Lagrangian Lq 
in l(32|) is equal to 2. This contradiction proves the statement. 

Finally we consider the example of a dynamical system for which an integrating multiplier and, consequently, 
the possibility of the Lagrangian description does not exist. Duglas [3] showed that the set of second-order 
equations 

x + y = 0, 
y + y = 

does not admit an integrating multiplier. Let us prove this. To this end, let us assume the opposite, namely, 
let there be a non-degenerate matrix hij that obeys equations ([9]), ifTSi) . (fT7|) . ifTS)) and l(23|) . Then from the 
algebraic equation (|23l) it follows that hij must be diagonal (/ii2 = /121 = 0), since in this case 

2 

Then, from condition (fT5|) we immediately obtain hn = 0, and arrive at a contradiction, det hij = 0. 

Thus, we can see that an action functional in second-order formalism cannot be constructed for some sets 
of differential equations. Nevertheless, as we show in the following section, it is always possible to construct an 
action principle for the equivalent set of first-order equations. 



3 Action principle in the first-order form 

Let a system with n degrees of freedom be described by a set of n non-Lagrangian second-order differential 
equations of motion. To construct an action principle, we replace these equations (which is always possible 
to do by introducing n additional variables, e.g. pi = if ) by an equivalent set of 2n first-order differential 
equations, solvable with respect to time derivatives. Suppose such a set has the form 

= r{t,x), (39) 
= {q\pi) , a = 1, ..,2n, 

where /" (t, a;) are some functions of the indicated arguments and by dots above we denote time derivatives of 
coordinates. Let us construct an action principle for this set. If (|39|) cannot be directly identified with Euler- 
Lagrange equations, then one can find an integrating multiplier, i.e. a nonsingular matrix which reduces the 
initial set of differential equations l(39ll to a variational derivative: 



<7„ [t] = {x^ - fit, x))^j^=0. (40) 



Since ga [t] is a variational derivative it must obey the Helmholtz criterion [T] 

Sga [t] ^ Sgf) [s] 
5x0 (s) (Jx" (t) ' 



(41) 



We will use this condition to find an integrating multiplier fi. In the general case one can assume that 
depends on time t, coordinates x" and time derivatives up to order to (m is a natural number), i.e. ga [t] = 



^We denote by Q the integrating multiplier for the first-order equations. 
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ga {t,x, ...,a;''"'). Having in mind this form of ga one rewrites l(4T|) as 

3=0 



5xP{s 

Differentiating the identity f {t) S (t ~ s) — f (s) S {s ~ t) over t one finds 

k 



/ (.) 5^=) (. - 1) = (-1)'= Cif^''^ it) '5^'"'^ - ^) ' = TlAv • (43) 



i=0 

Using l|43p we rewrite (|4T|) as 

4=0 j=Q 1=0 



dV dgi3 [t] 



SU~i) (t-s) . (44) 



Comparing in (|44l) the coefficient for (J^*^^ (i — s) one gets equations for g^ (t, x, ...,x^™^). When A: = we have 

^_V(_iyY-V^^ = o (45) 



Since ga [i] depends only on derivatives up to order m, the coefficient of the higher derivative ^"(^m) equation 
gSl) must vanish. This coefficient is (-1)"' 92g^/a2;"(™)aa;7(™) ^ which means that ga [t] must be linear on the 
derivatives of order m, i.e. 



go. \t\ = aa0 (t, X, X^"-!)) + ba {t, X, ...,X^^~^^) 



where Oq/j and ba are some functions. Since is a nonsingular matrix, (|40l) should be a system of first-order 
equations, i.e. we have m = 1 and = (t, a;). 

Now comparing the coefficient of 5^^'' {t — s) in l(44|) . we obtain: 

^af3 = -^I3a ■ (46) 

Then from l(45|) we have 

(9/3 (ila^/'^) - 5a {Vlp^p) + dt^aP + i'^ (S/jfia^ - da^p^i + 97^^/30) = . 

Since does not depend on x one gets the following equations for fi: 

da^f3'y + S/jfi^Q + d-y^ap = (47a) 

and 

dtnat} + If^afi^O, (48) 

where £ fflafs is the Lie derivative of fia/j along the vector field /'''and da = dj dx" , dt = d/dt. Thus we see 
that for a set of first order equations an integrating multiplier is a nonsingular matrix which depends only on 
time t and coordinates x", and obeys the conditions l(46|) . I|47ap and l(48l) . 

Let us analyze equations (|46l) - l(48|) for the matrix ilap following our work It is known that the general 
solution Q.ai3 of equation l(48|) can be constructed with the help of a solution of the Cauchy problem for equations 
(|39l) . Suppose that such a solution is known: 

= (i , a;(o) ) , xfo) = V'" (0, X(o) ) (49) 
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is a solution of equations l(39|) for any a;(o) = (^(o)) ' ^^'^ x"(^i^) is the inverse function with respect to 
'/f"(i,a;(o)), i.e., 

Then 

no.pit,x)^d^x''i^%Hx) dpx\ (51) 



where the matrix fl^^^ is the initial condition for flap, 



One can see that by choosing the matrix fl^^^{x) to be nonsingular and subject to the Jacobi identity, we 
guarantee the fulfilment of the same conditions for the complete matrix flafi (t, x) , since components of the 
latter are given by a change of variables fSTj) . 

Thus, we see that for any set of first order equations l(39|) . an integrating multiplier always exists, i.e. there 
always exists a Lagrangian L{t, x,x) which has the set of equations 

nap{x^ - f{t,x)) =0 (52) 

as its Euler-Lagrange equations. It is easy to see that this Lagrangian is Hnear in the first-order derivative 
since equations l(52l) do not contain second-order derivatives, i.e. the corresponding term d'^L/dx°'dx^ vanishes. 
The general form of this Lagrangian is 

L = J^i" - H , (53) 

where Jq — Ja{t,x) and H — H{t,x) are some functions of the indicated arguments. The Euler-Lagrange 
equations corresponding to l(53|) are 

^ = ||-^|| = 0=^ -daH - dtJa + {daJp - dpJa) x^ = 0. (54) 
Comparing equations l(52|) and fM)) one getfl 

^ daJp - dj}Ja , (55) 

and 

n^pf -dtJa^daH. (56) 

The functions Jq and H can be found from conditions l(55|) and l(56|) if the matrix fla/3 is given. One can see that 
consistency conditions for these equations exactly give us equations (|46l) - l|48p for an integrating multiplier flap- 
We recall that the general solution Ja{t,x) of the equation l(55|) . provided that flaf3 is a given antisymmetric 
matrix that obeys the .Jacobi identity, is given by 



Ja{t,x) = I x'^fli3ait,Sx) sds + da(pit,x) , (57) 

Jo 

where ip{x) is an arbitrary function. Substituting fSTj) into l(57|) . we obtain 



Equation l(58|) describes all the ambiguity (arbitrary symplectic matrix fl^^j and arbitrary function (^(t,x)) in 
constructing the term Ja{i, x) of the Lagrange function l(53|) . 

To restore the term H in the Lagrange function (|53|) , we need to solve the equation (jSG)) with respect to H. 
To this end, we recall that the general solution of the equation dif = gi, provided a vector gi is a gradient, is 
given by 

f{x) = I ds x'-gi{sx) + c, 
Jo 



^As a more abstract argument in favour of such a form of f! we invoke the Poincare Lemma, according to which any closed form 
is locally exact, but here we do not treate the global problems wich can arise from the non-trivial topology of x"-space. 
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where c is a constant. Taking the above into account, we obtain for H the following representation: 



1 



H{t,x)^ dsx^ [npa{t,sx)f°'{t,sx) - dtJp{t,sx)]+c{t) , (59) 







where c{t) is an arbitrary function of time, and flpa and Jg are given by (|5T1) and (|58| respectively. All the 
arbitrariness in constructing H is thus due to the arbitrary symplectic matrix fl'^^g and the arbitrary functions 
ip{t,x) entering the expressions for il^a and Jp, and the arbitrary functions c{t). 

We can see that there exists a family of Lagrangians l(53|) which lead to the same equations of motion l(39|) . 



It is easy to see that actions with the same ri^^°^ but different functions ip{t, x) and c{t) differ by a total time 
derivative (we call such a difference trivial). A difference in Lagrange functions related to different choice of 
symplectic matrices flf^ is not trivial. The corresponding Lagrangians are known as s-equivalent Lagrangians. 
As an example, let us consider a theory with equations of motion of the forrrH 

x^ A{t)x+j{t). (60) 

An action principle for such a theory can be constructed (see [15]) following the above described manner. 
The solution of the Cauchy problem for the equations l(60|) reads 

x(i) =r(t)x(o) +7(i), (61) 

where the matrix r{t) is the fundamental solution of ([60l) . i.e., 

t^AT, r(0) = 1, (62) 

and j{t) is a partial solution of ((6Ql) . Then following ijSTj). we construct the matrix 

n = A^n^°'>A, A = r-^ (63) 

and find the functions J and H according to (|58l) and ((59|) , 



J = ^xD. , H ^ ^xBx - Cx , (64) 

where 

B = ^{nA- A^n) , c = nj . (65) 

Thus, the action functional for the general quadratic theory is 

S[x] ^\ j dt {x^lx - xBx - 2Cx) . (66) 

In conclusion we note that it is always possible to construct a Lagrangian action for any set of non-Lagrangian 
equations in an extended configuration space following a simple idea first proposed by Bateman [22] • Such 
a Lagrangian has the form of a sum of initial equations of motion being multiplied by the corresponding 
Lagrangian multipliers and new variables. The Euler-Lagrange equations for such an action contain besides the 
initial equations some new equations of motion for the Lagrange multipliers. In such an approach one has to 
think how to interpret the new variables already on the classical level. Additional difficulties (indefinite metric) 
can appear in course of the quantization |23)-|26|. 

Acknowledgement 1 Gitman is grateful to the Brazilian foundations FAPESP and CNPq for permanent 
support; Kupriyanov thanks FAPESP for support. 



^Here we use matrix notation, x = {x°') , A{t) = (A(t)J^) , j{t) = (jit)") , a,f3 = 1, ..,2n. 
^For simplicity we choose the matrix n''^' to be a constant nonsingular antisymmetric matrix. 
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